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In the perturbation theory, trasnport phenomena in hot 
non-Abelian gauge theories like QCD are often plagued with 
infrared singularities or nonperturbative effects. We show, in 
the context of the Kadanoff & Baym formalism, that there 
are certain nonequilibrium processes which are free from such 
difficulties. For these processes, due to an interplay between 
the macroscopic and microscopic physics, characteristic time 
scale (the mesoscale) naturally enters as an infrared cutoff and 
purely perturbative description by the Boltzmann equation is 
valid. 
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Recently, there have been considerable interests in the 
transport phenomena in hot gauge theories moti- 
vated by the physics of the quark-gluon plasma (QGP) 
[ 3 1 and the baryon number violation in the early universe 
\7\. The standard approach is to derive the Boltzmann 
equation and obtain various transport coefficients from 
it. However, in this direction one often encounters log- 
arithmic singularities due to the exchange of very soft 
gluons. Such singularities are typically summarized in 
the following integral 



dq-, 
Q 



(1) 



where g is the gauge coupling constant and T is the tem- 
perature. The upper limit gT is needed to be consis- 
tent with the hard thermal loop approximation || lead- 
ing to the integrand 1/q. Due to this logarithmic sin- 
gularity, some quantities associated with nonequilibrium 
phenomena, e.g., the damping rate of thermal gluons and 
fermions, the color relaxation rate |22|,^3| are divergent. 
Henceforth we call these quantites transport coefficients, 
although it is an abuse of the nomenclature. In fact, they 
are superficially divergent. This is because the integrand 
1/q is valid only when q is larger than some nonpertur- 
bative scale \x. So we must cut off the integral at this 
scale anyway and further analyses require nonperturba- 
tive considerations, fi is considered to be of order g 2 T 
up to a possible logarithmic factor of g Q . 

The above situation seems to indicate that even in the 
lowest order of perturbation theory one cannot be free of 
the nonperturbative effects. In this paper we will show 
that this is not always the case at least in some physi- 
cally important processes like the damping of thermally 
excited particles. 

There are several derivations of the Boltzmann equa- 
tion. Here we follow one by Blaizot & Iancu jp based 
on the Kadanoff & Baym formalism H because, as we 



will see, it is well suited for our purpose. In this for- 
malism, the Boltzmann equation for the pure SU(-/V C ) 
plasma describes the slow variation of the Wigner trans- 
form G(k,X — ^5^) of the hard (containing spatial mo- 
menta of order T) gluon two point function G{x, y) 



2[k ■ D x ,SG(k, X)} - 2gk a F al3 {X)d f lG eq {k) 
= -T{k)5G{k,X) 



SE<G>, 



r(fc) = s< (k) 



(2) 



(3) 



where 5G = SG > = 6G < and 5E < are the deviations 
from equilibrium of the covariantized gluon propagator 
|l| and the Keldysh components of the gluon self en- 
ergy transversely projected, respectively. The covariant 
derivative D 1 ^ and the field strength F^ U {X) are con- 
structed from the supersoft (containing spatial momenta 
of order g 2 T ) mean field A»(X). 

Let us consider the problem of the damping of a hard 
gluon using (Q). To simplify (|J) one usually employs 
the Kadanoff & Baym ansatz which assumes that the 
nonequilibrium Green's functions are proportional to the 
spectral function. Furthermore, one assumes 

G > (k,X)~G < (k,X) =p{k.X) ~ 2ire(k )S(k 2 ), (4) 

where e(fco) is the sign function. This is known as the 
quasiparticle approximation. (Actually, this approxima- 
tion is required because when we have written down (j^) 
we dropped terms which produce the width in the spec- 
tral function S. Note that we have already utilized the 



fact that <5G ; 



5G < 



which is a consequence 



of the quasiparticle approximation.) Thus, one ends up 
with the following structure: 

5G{k,X) ee p(k, X)8N{k, X) 
~ 2TrS(k 2 )[d(k )6N{k, X) + 9(-k )SN(-k, X)}. (5) 

Dropping the mean field and neglecting the deviations 
from the equilibrium value other than the mode with 

wavevector k, that is, neglecting (5S < , we obtain 



dX, 



SN(k,X 



-2 7 <W(k,X ), 



(6) 



where 7 ee ^r(fco = k) is the damping rate of the hard 
gluon. It is well known that this quantity is divergent 
because of the long range nature of the color magnetic 
force. To see the divergence, we calculate it with 4- vector 
fc M slightly off-shell to obtain 



1 



T(k) ~4a£;T7V c ln 



feo — k\ 



(7) 



where u> p is the plasma frequency and a = |-. We see 
that the product of T(fc) and the on-shell delta function 
(|J) appearing in (|J) is clearly ill-defined. 

To avoid such a singularity, the first thing one can 
think of is to include the width to the spectral function. 
Perturbatively, the delta function approximation is quite 
a good approximation for a hard gluon. Nevertheless we 
try to include loop corrections. To the one loop order, we 
get the plasmon pole. T(k) is logarithmically divergent at 
the pole in the same manner as the bare case. In addition, 
we don't know how to handle with the Landau damping 
spectrum. To obtain the width, we must proceed to the 
next order. But the width is essentially the same as the 
damping rate which is divergent. So this approach is not 
practical besides the technical complexity. 

Another approach is to consider the finite time effect 
0. (See also (|||.) Notice that the delta function ap- 
pears as the Fourier transform of the propagators: The 
free and equilibrium propagators read 

G>(t,t') 

= 2-{(l + iV(fc))e- lfc ( t -*') - (1 + N{-k))e lk(t - tr) }, 
2k 

G<(t,t') = -L{jV(fc) e - 4fc (*-0 - N{-k)e lk ^^}, (8) 

where N(k) is the Bose distribution function. Fourier 
transforming with respect to the relative time from minus 
infinity to plus infinity, we obtain 

/oo 
d{t-t')e iko ^-^G^(t-t') 
-oo 

= 27re(k )S(k 2 )N(k ). (9) 

Now we perform a modification needed in the presence of 
infrared singularities. Consider a situation that at time 
t = we add a hard gluon with momenta k^ to the 
initially equilibrium plasma. We write the deviation from 
the equilibrium at t, t' > as (See (||).) 

6G<(t,t>) = SG>(t,t') = ^SNCk 1 X )e-^ t - t '\ (10) 

Once we have fixed Xq, t and t 1 are constrained as follows 



t, t' > 0, X = 



t + t' 



2X > t - t' > -2X . (11) 



Our point is that we perform the "Wigner transforma- 
tion" in this finite time interval 



6G(k,X ) 



2X„ 



d(t - t') e *(fco-fc)(t-0 



SN(k,X ) 



~2X Q 2fc 
5N(k 7 X ) 2sin[2X (fc - k)] 



2k 



kn — k 



(12) 



Thus the "spectral function" acquires the width ~ 1/ Xq 
by the finite time effect. Note that because SG^ = SG^, 
3) is compatible with the quasiparticle approximation 
Substituting ([l^) in (||), we get after fco-integration 



2k— SN(k, X ) -AakTN c \n(2uj p X )SN(k,Xo), (13) 
dX 



6N(k, Xq) ~ 6N(k, 0) exp{-2aTN c X a ln(2cj p X )}. (14) 

We see that the gluon distribution function exhibits the 
anomalous (non-exponential) damping 1213[|. The in- 
frared singularity is cut off by the inverse of the time 
scale with which we are looking at the system. Our 
method leading to (|l4|) is different from ||] in that we 
have included the finite time effect entirely into the spec- 
tral function. This is very convenient in the later argu- 
ment. 



It is clear that (|14|) is valid only for 1/gT < X < l/V 
This is because (0)is only valid when gT \k — k\ > 
\x and in the integral leading to (|l3|) the region Jfc — 
k\ ~ 1/Xq is important. However, in this region ( |l4| ) is 
a definite prediction of the perturbation theory and no 
nonperturbative effect comes into play. 

What we have done above can be stated in more gen- 
eral terms. Consider a nonequilibrium process (not re- 
stricted to gauge theories) which has a characteristic (or 
macroscopic) scale of the temporal variation. Charac- 
teristic scales of the spatial variation are less important 
compared to the time scale in the following argument. 
(Note that transport processes occur even in spatially ho- 
mogeneous systems.) We introduce a time scale r such 
that t is much smaller than the characteristic time scale 
and much larger than 1 / gT. Imagine a space-time lattice 
whose lattice spacings are given by r and characteristic 
spatial scales. (See, e.g., 0,|5).) Each cell is specified 
by the Wigner coodinate X^. We can "Fourier trans- 
form" only within a cell to discuss the difference between 
neighboring cells, that is, to discuss transport phenom- 
ena. Usually, scales of microscopic physics are literally 
much smaller than macroscopic scales so the integration 
ranges are practically infinite. But in the present case 
they cannot be sent to infinity. This effect practically 
smears the spectral functions making otherwise divergent 
transport coefficients logarithmically dependent on r. It 
is interesting to see that the infrared singularity originat- 
ing from the spatially long range color magnetic force is 
cut off by the time scale. Here we can see an interplay 
between the macroscopic and the microscopic physics, or 
the feedback effect of the macroscopic physics. Namely, 
the mass singularity coming from the microscopic colli- 
sion is remedied by the nonequilibrim nature of the bulk 
system whose evolution in turn is generated by the accu- 
mulation of microscopic collisions. 

In (Q), k is an arbitrary 4- vector, which enables us 
to study around the mass shell. Due to the finite time 
effect particles participating in a collision have certain 
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off-shellness. This effect can be naturally incorporated 
in the Kadanoff & Baym formalism |l6) . 

The full Boltzmann equation (0) becomes JD 

(v ■ D x ) ab W b {X,v) = v - E a {X) - j{\V a {X,v) 



dCl' 



(vv') 2 



} w a (xy) , (15) 



ttJ 4tt ^1 - (v-v') : 



where = (l,k) and SN(k,X) = -gW{X,w)^§. 7 is 
commonly called the damping rate of the hard gluon, but 
the logarithmic singularity is cut off by r, the intermedi- 
ate time scale between the macroscopic and microscopic 
physics. If we are interested in the short time evolu- 
tion of processes with the characteristic time less than or 
equal to 1/fi, we can safely use (|l5|) which is free of any 
nonperturbative effects ]17| ]. 

Note that the ordinary Boltzmann equation for a dilute 
classical gas can be used only for an intermediate time 
scale (the mesoscale) much longer than the collision time 
and much shorter than the mean free time as was first 
pointed out by Bogoliubov The situation is analo- 
gous here. turns out to be of the same order as the 
mean free time of hard gluons. So the purely perturbative 
Boltzmann equation proposed here is in accordance with 
Bogoliubov's spirit and describes intermediate processes 
leading to the damping of the hard gluon and the color 
transportation p9] . This is our main conclusion. 

The above consideration cannnot be extended to the 
situations when the characteristic time scale is larger 
than 1/fJ.. Beyond this regime, we must cut off the inte- 
gral by the momentum cutoff /j, instead of \ko — k\. So 
we must regard the resulting equation as an effective the- 
ory and further analysis requires some nonperturbative 
methods. One of the most interesting examples of this is 
the baryon number violation in the early universe. It was 
shown that the characteristic time scale of that pro- 
cess is l/g 4 T much larger than The effective theory 
was constructed by Bodeker Q and the hot sphaleron 
rate has been calculated by Moore on the lattice plfl . 

In summary, we have shown, in the Kadanoff & Baym 
formalism, that the nonequilibrium non-Abelian plasma 
may exhibit mesoscopic physics which can be described 
by the purely perturbative Boltzmann equation. Al- 
though we do not know phenomcnologically interesting 
processes that have such short characteristic time scales 
at present, these may be found, for example, in the ini- 
tial stage of the QGP in future ultrarelativistic heavy ion 
experiments. 
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